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Abstract 

We use information about DIS and J /ip production on hydrogen to model the f-dependence of 
the 7*iV scattering amplitude. We investigate the profile function for elastic scattering of hadronic 
components of the virtual photon off both a nucleon and heavy nuclear target, and we estimate the 
value of the impact parameter where the black body limit is reached. We also estimate the fraction 
of the cross section that is due to hadronic configurations in the virtual photon wave function that 
approach the unitarity limit. We extract, from these considerations, approximate lower limits on 
the values of x where the leading twist approximation in DIS is violated. We observe that the 
black body limit may be approached within HERA kinematics with Q 2 equal to a few GeV 2 and 
x ~ 10 -4 . Comparisons are made with earlier predictions by Munier et al., and the longitudinal 
structure function is compared with preliminary HERA data. The principle advantage of our 
method is that we do not rely solely on the t-dependence of /j-meson production data. This allows 
us to extend our analysis down to very small impact parameters and dipole sizes. Finally, we 
perform a similar calculation with a 208 Pb target, and we demonstrate that the black body limit 
is already approached at Q 2 ~ 20 GeV 2 and x ~ 10 -4 . 
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I. INTRODUCTION 



One of the current theoretical challenges in quantum chromodynamics (QCD) is to de- 
scribe high energy interactions with hadrons in terms of fundamental field theory. It is 
observed that high-energy hadron-hadron scattering interactions become completely absorp- 
tive (black) at small impact parameters so that elastic scattering can be viewed essentially 
as a shadow of the inelastic cross section in the sense of Babinet's principle. If this regime 
occurs at most of the impact parameters which contribute to the inelastic cross section, then 
the elastic and inelastic cross sections are equal. This limit is often referred to as the black 
body limit (BBL) in analogy with the quantum mechanical situation of scattering from an 
absorptive share of radius r in which case the total cross section is equal to 2nr 2 , (see e.g., 
problem 1 of section 131 in Ref . Q]). This limit is also loosely referred to as the unitarity 
limit, although unitarity alone admits cross sections as large as Airr 2 provided there are no 
inelastic interactions (see e.g., problem 2 of section 132 in Ref. [l|). I n this paper we will 
assume, in line with experiment, that the amplitude is predominantly imaginary. In this 
case, the unitarity limit coincides with the BBL. The black body limit in deep inelastic 
scattering (DIS) is an interesting new regime in QCD where the coupling strength is small, 
but where the leading twist approximation breaks down and new small coupling methods 
are needed. In order to understand the transition to the BBL, it is important to construct 
models which accurately describe 7* N scattering over a wide range of kinematic variables. In 
particular, such a model should interpolate smoothly between the nonperturbative domain 
and the region where interactions are accurately described by leading twist perturbative 
QCD (pQCD). In a recent paper by Munier et al. an estimate of the proximity to the 
BBL for the interaction of a color dipole with a proton was made using data from diffractive 
electroproduction of p-mesons 3- The techniques used in 3 are limited by the need to 
model the p-meson wave function, as well as by limited information on the t-dependence 
for p production. Within the impact parameter representation, this means that their pre- 
dictions are constrained to intermediate values of impact parameter (b > 0.3 fm) |3j. (It 
will be important to keep in mind the distinction between impact parameter, b, and dipole 
size, d.) We will find that our analysis is valid down to very small values of b and d. We 
will make an improved estimate of the onset of the BBL by combining information from 
DIS and J/ip production. A great advantage of our technique is that, unlike models which 
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are restricted to using meson production data, our model is valid down to very small dipole 
sizes because we utilize leading twist pQCD for calculations involving small size qq pairs. 
Within our model of 7*iV interactions, we always assume that the virtual photon can be 
written as a linear combination of hadronic states. Furthermore, we relate the hardness 
of the interaction to the size of the hadronic state as is done in Ref. 0], and we take into 
account the dependence of the hadronic interaction on the virtuality of the original photon. 
Large size configurations constitute the soft component of the interaction whereas small size 
qq pairs constitute the hard component. 

Modeling the 7*iV interaction requires making three principle observations. The first is 
to recall that the total cross section for small size dipole configurations {d < .3 fm) is known 
within the framework of perturbation theory down to x ~ 10~ 4 . Therefore, we will have a 
complete picture of the interaction of small size configurations with nucleons if we can extract 
the t-dependence from experimental data. In Sect. [TTl we discuss how the t-dependence of the 
small size qq pairs can be extracted from J /if) production data. Secondly, we note that the 
soft scattering of large size hadronic configurations is well understood phenomenologically in 
terms of effective Pomeron exchange. Hence, our model should reproduce the pion-nucleon 
amplitude for hadronic configurations comparable to the size of the pion (d ~ .6 fm). Finally, 
we must model the behavior of the amplitude for intermediate hadronic sizes (.3 fm < d < .6 
fm). This is an interesting and poorly understood region of kinematics, and our model will 
allow for readjustments in the transition region. 

In Ref. 0| , a system was devised for relating the transverse size of a qq dipole to the virtu- 
ality, Q 2 . The pQCD result for the inelastic 7*iV cross section was interpolated to large size 
hadronic configurations and matched to the total cross section for pion-nucleon scattering. 
The hadronic size, d, is used to interpolate between the hard and soft regions in this paper 
as it was in Ref. Jj]. d represents the transverse size of a quark- ant iquark pair in the limit 
that d is very small {d < 0.1 fm). However, as d grows large, the dipole picture becomes 
inappropriate since the hadronic components that correspond to soft interactions consist of 
large, complex hadronic states. Larger values of d should be interpreted as transverse sizes 
of general hadronic components of the virtual photon wave function. The cross section takes 
the form familiar from perturbation theory: 




(1) 
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Here, A = Q 2 d 2 is a universal scaling ansatz used to relate energy scales, Q 2 , to trans- 
verse dipole sizes, d; x is the light-cone fraction for the gluon attached to the qq loop. 
One manifestation of the QCD factorization theorem is that the contribution of hadronic 
configurations within the photon to the longitudinal cross section, <7l, is peaked around a 
narrow range of small dipole sizes (see Fig. 0J). The value of A is chosen so that, within the 
perturbative region, d 2 = X/Q 2 is approximately the average dipole size contributing to ol- 
For large Q 2 , A takes on values of the order of 10. In fact, it is found that F 2 and Fl depend 
very weakly on the value of A within the perturbative region of d 4]. Changing the value of 
A thereby provides a universal parameter for tuning the cross section within the transition 
region. In this paper, we use A = 4 because it is found that A = 4 best describes J /if) 
data over a wide range of kinematics (both perturbative and nonperturbative) . Figure El 
compares the profile function for A = 4 and A = 10. The function used in Ref. J] to inter- 
polate between the hard and soft regions matches smoothly to the pQCD result at small d 
and to the pion-proton cross section at large d. Furthermore, it takes into consideration the 
breakdown of the leading twist formula in the small x limit. 

Notice that the qq pair is not a fundamental object since it is always produced off-shell 
by a virtual photon. This is taken into account on the right side of Eq. through implicit 
dependence of — , on Q 2 and Q 2 (see Ref. |4j and Eq. The dependence of the dipole 

cross section on the external photon virtuality, Q 2 , is a feature that is absent in other models 
such as the one proposed by Wiisthoff and Golec-Biernat ,5]. This point will be important 
for what we discuss later because it means that we cannot speak unambiguously about the 
dipole cross section with referring to the interaction for which it is a sub-proccess. 

In this paper, we will model the t-dependence of the 7*iV elastic scattering amplitude 
(which, of course, cannot be observed experimentally) using data from J /if) photo(electro)- 
production in conjunction with the pion-proton elastic scattering amplitude. As when we 
model the behavior of the total cross section in various kinematic regions, we model the 
t-dependence of the j*N amplitude by considering three distinct steps. First, we model the 
t-dependence of the small size qq configurations. Because of transverse squeezing in the J/ if) 
wave function, data taken from J /if) production is appropriate for use in modeling the in- 
dependence of small size, hard scale qq configurations. Next, the t-dependence of soft, large 
size hadronic configurations is approximated by the pion structure function, the pion being 
a reasonable approximation to a large size hadronic component of the virtual photon wave 
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function. For soft interactions, a factor for soft Pomeron exchange is included to account 
for the slow rise in cross section at small x. Finally, we use the transverse size, d, as a 
parameter to interpolate between hard and soft physics. In our analysis, we will transform 
our expression for the amplitude into the impact parameter picture to look for the regions 
where the impact parameter space amplitude approaches the unitarity limit, and thus to 
estimate the values of the impact parameter where the BBL is attained. 

It should be noted that our model has limitations which restrict how it can be applied. 
It is important to keep in mind that Eq. ((TJ is multiplied by a color factor of 9/4 when 
the hadronic configuration is a color singlet composed of octet representations of SU(3), as 
in a gluon dipole. Such interactions are expected to be abundant in the small x limit, so 
the BBL will be reached at larger impact parameters than what is predicted by considering 
only the interactions of small size qq pairs. Furthermore, at very small x, effects from 
cross section fluctuations of the virtual photon become important, and taking into account 
only elastic dipole scattering becomes inappropriate. In particular, the total cross section 
in hadron-hadron scattering has a significant contribution from inelastic diffraction. In 
Ref. Miettinen and Pumplin write the contribution from inelastic diffraction in terms 
of fluctuations around hadronic eigenstates. Estimates of the contribution from hadronic 
fluctuation, to the cross section snow that it is uot negligible (see, e.g. Ref. Q). lu fact, with 
a decrease of the dipole size, the relative importance of inelastic diffraction increases. The 
importance of inelastic diffraction for small qq sizes is discussed, for example, in Ref. Q. The 
important point to note here is that inelastic diffraction will contribute to the breakdown 
of the leading twist approximation at low x before the BBL for elastic scattering of the 
hadronic configurations in the photon wave function is reached. As such, we do not seek 
to place an absolute boundary on the region where corrections to the DGLAP evolution 
equation become relevant. Rather, we construct a model that puts a lower limit in impact 
parameter space on regions approaching the BBL. Furthermore, since the leading twist 
approximation is not accurate in the vicinity of the unitarity limit, the BBL establishes a 
lower limit in impact parameter space on regions where the DGLAP equation is applicable. 

In Sect. |nj we outline our model for the t-dependence of the hadronic configuration- 
nucleon amplitude. In Sect. Hill we transform our expression for the amplitude into the 
impact parameter representation and study the proximity of the profile function to the 
unitarity limit as a function of impact parameter, b. We write the transverse and longitudinal 
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cross sections, crj L , as convolutions of the dipole cross section with the photon wave function 
in Sect. IIVI Section IIVI concludes with a calculation of the fraction of the 7*iV cross section 
due to large values of the hadronic profile function. In Sect. EJ we compare our results to 
an earlier study of the S-matrix t-dependence in impact parameter space. In Sect. IVII we 
perform the same calculation for the situation where the target is a 208 Pb nucleus. For the 
case of a nuclear target, we match the pQCD calculation at small d to the Glauber multiple 
scattering theory result at large d. Finally, we summarize our observations in the conclusion. 



II. MODELING THE ^DEPENDENCE 

Starting with the expression for the total cross section in Eq. (JTJ), we devise a model for 
the scattering amplitude by writing it in the form, 

A hN {s, t) = isa tot f{s, t) , (2) 

where f(s,t) accounts for the t-dependence of the interaction, and atot is determined from 
the QCD improved dipole picture (Eq. (JTJ). The "hat" on a tot is to distinguish the total 
cross section for the scattering of one component of the photon wave function from the total 
7*iV cross section which we consider in Sect. IIVI Applying the optical theorem in the large 
s limit reproduces atot- For now we assume that the amplitude is purely imaginary. We will 
return to the question of a real part of the amplitude at the end of Sect. 11111 

In this section, we will make an estimate of the form of f(s, t) which will take into account 
the nonzero size of the qq dipole and which will smoothly interpolate between perturbative 
and nonperturbative regimes using the hadronic size, d, as a parameter. We act in the spirit 
of Ref. ^ by modeling the i-dependence of the amplitude in the soft and hard regimes and 
by using d to build a smooth interpolation. The three steps: building a model for the small 
dipole region, building a model for the large wave-packet region, and interpolating between 
the two regions are outlined in the next three paragraphs. 

We start by writing the general structure of the amplitude. The t-dependence, f(t, x, d), 
is written as the product of three functions 

f(t, x, d) = F N (t, d)F h (t, d)F P (t, x, d) . (3) 

Here and in the rest of this section, the dependence of / upon s is replaced by dependence 
upon x and d. F N (t,d) describes the t-dependence of the nucleon target, F h (t,d) describes 



the t-dependence of the hadronic projectile, and Fp(t,x,d) accounts for Gribov diffusion. 
This method of separating the t-dependence into three factors corresponding to different 
sources of t-dependence is similar to what is used in Ref. 

The next task is to model the small dipole size t-dependence. Both the soft Pomeron 
exchange factor and the hadronic form factor must approach unity as the size, d, shrinks 
to zero. The QCD factorization theorem implies that the t-dependence of the small dipole- 
nucleon amplitude is universal. Hence, it can be extracted directly from J/ip photo[electro)- 
production since the J/ip wave function is known to be a small size wave-packet [lfj. Data 
from J/ip production reveal that the two-gluon form factor is, 

F„M-o) = F 1 ( ( )~ lr -i^. <«) 

The subscript, 1, labels the two-gluon form factor and m\ is a measurable parameter in the 



two-gluon form factor. The value, m\ « 1.1 GeV 2 is extracted from data in Refs. |llLll2l. 13J. 
For a detailed discussion of the two-gluon form factor in J/ ip production, see Ref. {14] . In 14] , 
it was discussed in detail how, due to the transverse squeezing of the J/ip wave function, the 
t-dependence of the amplitude comes solely from the two-gluon form factor. In particular, 
it was found that the dipole form factor contributes only about 0.3 GeV~ 2 to the slope 
of the t-dependence. The assumption that only the gluon form factor is relevant for J/ip 



production has been successfully tested against data in Refs. (2Jj Ilia, ll2( . Hence, in the limit 
of small dipole sizes, 

f(t,x,d^0)=F l (t). (5) 

Next we construct a model for the large wave packet behavior. When the hadronic 
state has a large size, the t-dependence receives contributions from sources other than the 
two-gluon form factor. We rewrite Eq. Q in the form, 

/(*, x, d) = F e N m -(t, d)F h (t, d)F P (t, x, d). (6) 

Now, Ffj m {t) is the electromagnetic form factor of the nucleon which is known phenomeno- 
logically to take the form, 
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where Wq pa 0.7 GeV 2 . Large size hadronic configurations can be reasonably expected to 
have t-dependence similar to the pion electromagnetic form factor. Thus, for the hadronic 
form factor we use the well-known form of the pion form factor, 

with m\ pa 0.6 GeV 2 . Here, d n is the characteristic size of the pion, and takes on a value of 
approximately 0.65 fm. This value for the pion size is consistent with what is used in the 
matching ansatz of Ref . J^j and agrees well with data for the 7rN cross section in Ref . [lij . 
For low-x soft scattering there is also a factor that arises from Gribov diffusion effects: 

F P (t,x,d^ d n ) ~ e - a ' tln ^ . (9) 

The factor, F P (t, x, d), describes the exchange of a soft Pomeron with Regge slope a' pa 0.25 
GeV -2 and Xq = 0.01. The value of Xq is determined by the boundary of the region where 
Gribov diffusion effects become significant. 

Finally, we must find a reasonable way to interpolate between the hard and the soft 
regions. We use the t-dependence discussed in the previous two paragraphs to guess the 
following form for the hadronic configuration-nucleon amplitude: 



A hN (s,t) = isa tot {l _ t/ l M2{d2))2 1 _ td l /d> f • (10) 

In order to give the variation with d geometric behavior, we use d 2 as a parameter. To in- 
terpolate between the nucleon and the two-gluon form factors, we have defined the function, 

MV)H mS " K " mS) ^' d - d ' ■ (ID 

777.Q , otherwise 

Note that when d equals d n , AhN is the product of Eqs. (JZJ), (jHJ)) and (JUJ). In the small d 
limit, the dipole form factor and the Pomeron form factor approach unity, M 2 (d 2 ) — > raf, 
and the limit in Eq. (jHJ) is recovered. Varying d 2 interpolates smoothly between the hard 
and soft regions. Note that we neglect a possible small x dependence of F N (t, d) at x < 0.01. 
(See the discussion in Ref. [19].) However, our model is adjusted to reproduce the observed 
x-dependence of the slope for photoproduction of J/ip mesons. 
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III. IMPACT PARAMETER ANALYSIS 



Having obtained Eq. (fTU|) . the next step is to transform to the impact parameter repre- 
sentation where the profile function is defined by the relation, 

A dipole , N (s, t) = lis J d 2 be-^ rb T h (s, b) , (12) 

where t = —q 2 . The subscript, h, indicates that we are considering the profile function for 
the scattering of a single hadronic component of the photon wave function from the proton. 
We get the profile function by inverting Eq. (fL2"|) . 

r h (s,b) = — L- J d 2 qe^ b A hN (s,t) . (13) 

For an imaginary amplitude, the BBL is reached when I\(s,fe) = 1 and the elastic and 
inelastic cross sections are equal. Recall that if a singlet dipole consists of color octet 
representations of SU(3), Eq. has an extra factor of 9/4, so the BBL for the interaction 
of a hadronic configuration with the nucleon is certainly reached for Th{b) ~ 1/2. Therefore, 
whenever > 1/2, we conclude that the interaction takes place near the BBL. In the rest 
of this section, we will suppress explicit reference to the argument, s, in the profile function. 
We have plotted the function T/ l (6) for different values of the dipole size and x in Fie 
We have used gluon distributions from CTEQ5L in the perturbative calculation of &tot |21 
Recall from the introduction that our model requires that we specify the external photon 
virtuality. Since we are interested in the possibility of reaching the BBL at a few GeV 2 , we 
have set Q 2 = 2 GeV 2 in Fig. □ 

A Gaussian ansatz is commonly used in experiments to extrapolate the t-dependence to 
large values. Let us, therefore, briefly compare the behavior of our model to that of a simple 
Gaussian. We start with the form 

A hN {s,t) = isa tot e ct/2 , (14) 

which is then transformed into impact parameter space giving, 

= |^ . (15) 

The slope of the Gaussian, C, is chosen so that it yields the same standard deviation in r^(6) 
as our model. One danger in using a Gaussian model is that it neglects the importance of 
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FIG. 1: The hadronic configuration- nucleoli profile function for different x values. The large 
region ( > 1/2) is reached for intermediate hadronic sizes. (See figure [3) Here, Q 2 is taken to 
be 2 GeV 2 . 

interactions in peripheral regions. Our model attempts to fix this problem by spreading out 
the distribution in t. Note in Fig. Elthat our model fall off more slowly with b. 

Now let us estimate the contribution of large values of T h (b) to the total hadronic cross 
section. The total cross section follows from the optical theorem, 

a tot = 2 J d 2 bReT(s,b). (16) 

We have made a numerical evaluation of the fraction of the total hadronic configuration- 
nucleon cross section obtained by setting different upper limits on the 6-integral in Eq. (|16|). 
In Fig. El one can see that no more that about thirty percent of the total hadronic cross 
section is due to values of Th > 1/2. Moreover, contributions from large values of 1^(6) 
occur for hadronic sizes close to the pion size, d ~ .6 fm. Averaging over the photon wave 
function will lead to a suppression of contributions from larger size hadronic configurations, 
so there will indeed be a small contribution to the total DIS cross section due to large values 
of T h {b) (see Fig. HJ). The goal of Sect. |TV|will be to determine whether the contribution to 
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FIG. 2: Comparison of the b behavior for our model with that of a Gaussian model. Our model 
falls off more slowly with b. The slope of the Gaussian used here is 0.17 fm 2 . 

the 7*iV cross section from large values of 1^(5) is significant enough that we may expect to 
see black body behavior within HERA kinematics. 

To summarize, Fig. ^ demonstrates that large values of Th(b) are approached for hadronic 
configuration-nucleon scattering at central impact parameters, b < 0.5 fm. In Fig. ^ it is 
seen that this is particularly true for hadronic sizes around d ~ 0.6 fm. Figure El shows that 
for d ~ 0.6 fm, a maximum of about 1/3 of the total hadron-nucleon cross section comes 
from values of Th(b) that approach the black limit. When d < 0.2 fm, a very small fraction 
of the total hadronic configuration-nucleon cross section comes from large values of Th{b). 
The only contribution from > 1/2 to the total cross section for d < .2 fm occurs at very 
small x (x < 10~ 4 ). 

Most of the model dependence in this calculation comes from uncertainty in the large-t 
behavior of the amplitude. The different curves in Fig.EJdemonstrate how our model changes 
if we remove contributions from large t. Notice that simply removing the contribution from 
—t > 3.3 GeV 2 leads to an error of less than ten percent. Thus, we do not expect our 
uncertainty in the large t behavior to have a drastic effect. 
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FIG. 3: The fraction of atot with contributions coming from values of greater than the corre- 
sponding point on the x-axis. 



We should also remark that we have considered only the non-spin-flip interactions. Cor- 
rections which account for the spin-flip amplitude would result in a smaller non-spin-flip 
amplitude than what we consider here. Experimental results in Ref. [2]] demonstrate that 
the polarization, P, is less than 0.2 for the range of t we are discussing. From the formula 
relating P to the spin-flip amplitude, 

2Im{A ++ A* + _) 2\A + _\ 



P 



I A, 



I A, 



I A 



(17) 



we find the fraction, !^~| < 0.1. Here, A ++ represents the amplitude with no spin-flip 

whereas A^ represents the amplitude with spin-flip. 

We now return to the issue of a real part of the amplitude which we ignored in Sect. 
In the considered kinematic region (t < —2 GeV 2 ), the ratio of the real to imaginary part 
of the amplitude, 7], is rather small. Indeed, if we adopt power law behavior for the total 
cross section, a to t ~ s p , we can estimate the value of 77(0) using the following formula which 
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FIG. 4: The distribution of the integrand in Eq. Q26|) over hadronic sizes for both the transverse 
and longitudinal cross sections. 

follows from the Gribov-Migdal result jz| at high energies in the near forward direction, 

ReA hN (s,t) 7id\na tot 



n(t) 

/w ImA hN (s,t) 
The amplitude can be rewritten as, 



7T 



2 dins 



A hN {s,t) -> s{i + r]{t))f{s,t), 



(18) 



(19) 



where the function, f(s,t), is assumed to be strictly real. For soft regions, the total cross 
section has the approximate s-behavior of the irN cross section as in Ref. Q, consistent 



with the behavior of a Donnachie-Landshoff soft Pomeron |23|]. In that case, p ~ 0.08, and 
Eq. (JT5j) gives r\ ~ 0.1. The second term in Eq. (fl"9"j) appears squared in the calculation of 
the cross section, so the correction to the cross section is approximately one percent. For 
the high Q 2 , low-x region, the total cross section experiences rapid growth and p ~ 0.25, or, 
by Eq. (|18jl. rj w 0.35. The correction to the squared amplitude is therefore approximately 
ten percent near the forward direction. Away from the forward direction, one must account 
for the small variation of rj with t. The effect can be estimated by considering the signature 
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FIG. 5: A demonstration of the rapid convergence of the profile function. Here, the profile function 
is plotted for different values of the upper limit, U, on the integral over t (— t = U 2 ). 

factor in the general form of the Reggeon amplitude. For —t < 2.0 GeV 2 , i](t) continues to 
contribute a negligible amount to the amplitude. 

The elastic cross section associated with Eq. (flTJJ) is found by integrating the profile 
function over impact parameters, 



a el = J d 2 b\T(s,b)\ 2 . (20) 
The total cross section is found using Eq. ()16|) . Therefore, the inelastic cross section is, 

Zinei = J d 2 b(2ReT(s,b) - |r( S ,6)| 2 ), (21) 
with the unitarity constraint, 

2ReT(s,b) - \T(s,b)\ 2 < 1. (22) 

If the amplitude is purely imaginary, then a e i < &i ne i and the unitarity constraint is that 
Fh < 1. Note that by considering only the imaginary part of the amplitude, we have 
considered only the real part of the profile function. If the amplitude is given a real part 
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FIG. 6: Comparison of the profile function for different values of A in the case of an intermediate 
hadronic size equal to 0.4 fm. Th changes by about fifteen percent at b = if A is changed from 4 



to 10. 



correction, then the profile function will obtain an imaginary part, the elastic cross section 
will increase, and the inelastic cross section will decrease. The correction to the unitarity 
constraint on ReT is, from Eq. (|22|). — (i]Re(T(s, b)) 2 . In the region of large Q 2 , the effect 
of a real part in the amplitude would clearly be noticeable. By Eq. ()22|). the unitarity limit 
on the real part of the profile function for t] « 0.35 would be, 



ReT < ~ 0.9. (23) 

~ 1 + T] 2 K ' 



Thus, the unitarity limit on the real part of the profile function may be less than unity 
by as much as ten percent. At this point, we should remark that both the contribution from 
the real part of the amplitude and the contribution from inelastic diffraction will tend to 
raise the boundary in impact parameter space where the BBL is reached. We neglect both 
effects in our model. As a consequence, when our model predicts that the BBL has been 
reached below a certain impact parameter, we can be confident that the same would be true 
in a model that incorporates inelastic diffraction and the effects of a real component of the 
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FIG. 7: These graphs are identical to those in Fig. ^ except that the value of Q 2 used to make 
each graph is calculated from the hadronic size. Note the larger values of the profile function at 
small d compared with Fig. ^ 

amplitude. On the other hand, if our model predicts that the BBL has not been reached, 
we must keep in mind that corrections due to inelastic diffraction and a real part of the 
amplitude may be important. In other words, the BBL may already be approached at larger 
values of b than what our model predicts. 

The uncertainty in the matching region is expressed by the uncertainty in the parameter, 
A. However, values of the order of 4 to 10 seem to work well and, as shown in Fig. there 
is a variation of only about fifteen percent at small impact parameters when we vary A from 
10 to 4. Note that this is done for a hadronic size of 0.4 fm which is in the region where 
the dependence upon A should be at its greatest. However, there remains another subtlety 
related to the matching of kinematic regions. First, recall the distinction between the energy 
scale, Q 2 , denoting the virtuality of the photon in a particular scattering process, and the 
scale, Q 2 which is the energy scale related to the hadronic size, d, through the scaling ansatz 



of Ref. [4j. These two scales are nearly equal as long as we consider hadronic sizes in the 
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FIG. 8: Demonstration of reasonable agreement between the color dipole model and recent HERA 
data 26] for F 2 at low Q 2 . The diferent_ curves correspond to the different parton distributions 
CTEQ6L,CTEQ5L, and MRST98 



2(j, 2^> 3- I" ° ,u ' calculations we used CTEQoL parton distri- 



butions because this yields optimal agreement between the dipole model and current data. 



vicinity of the average hadronic size for F 2 . In determining which value of x (called x in 
Eq. JJJ) should be used to calculate the hadronic cross section, the authors of Ref. 4] chose 
to relate the value of x to the value x for a particular 7*iV process in such a way that x 
varies as d~ 2 and so that for typical hadronic sizes, (x) = x . Within the color dipole picture, 
this accounts for the dependence of the the light-cone fraction, x , of the gluon attached to 
the qq pair on the mass, M 2 , of the dipole: 

Q 2 + M 2 



x 



(24) 



The result is that we cannot speak unambiguously of the hadronic cross section without 
referring to the virtuality of the probe which generated a component of given transverse 
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25]. 



size. From Ref. |4| we have, neglecting the constituent quark mass, 



O 2 

x =a;(l + 0.75gr). 

Q 



(25) 



Here we see that if we consider a fixed Q 2 , the universality of the hadronic cross section fails 
for small hadronic sizes (large Q 2 ), but is recovered for larger hadronic sizes. The value of 
x used in a calculation of the hadronic cross section will be significantly larger than x for 
small hadronic sizes, leading to a suppression of the cross section in the small size region. In 
particular, in Fig. the approach to the BBL at small d is slowed due to the large values of 
x needed to push the small size configuration on shell. In investigating the hadronic profile 
function, it may also be reasonable to determine Q 2 by letting it equal Q 2 so that the value 
of d always corresponds to a typical component of the virtual photon. We have done this 
in Fig. [7| and we can see that at small d, T h (b) is substantially larger, especially at small 
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x. Comparing Figs. ^ and [71 we see that at d = .1 fm this effect is significant while at 
intermediate hadronic sizes the effect is very small. For d > .5 fm there is no discernible 
difference between the two cases. The physical meaning of this effect is that the profile 
function for a small size configuration approaches the BBL more slowly if it is far off shell 
for a given Q 2 . Note that once we begin to calculate the total cross section, an external 
value for Q 2 is explicit, and we no longer have this ambiguity. 

Furthermore, there is some uncertainty in the gluon distribution used to calculate a to t- 
This is demonstrated in Figs. |H1 and where we compare results for the structure functions 
usine; CTEQ5L H, CTEQ6L Q, and MRST98 |2j| leading order gluon parton distri- 
butions. The dependence upon the parton distribution is seen to be small, but we used 
CTEQ5L parton distributions for all other calculations because they seem to yield optimal 
consistency with data. 

As we mentioned in Sect. UTl the value of gL that we used is consistent with the slope of 
the „N c rOSS S ect 10n as _d ln Re, Q and with the ma tch mg mS at Z US e d in Re, Q. 
In the model of the t-dependence, d n determines where soft Pomeron behavior becomes 
important, and one may well ask whether a different value of d n is appropriate. For models 
with a larger value of d n , the suppression of the profile function due to the Pomeron form 
factor, Fp, does not occur until one considers larger hadronic configurations. Therefore, 
for intermediate hadronic sizes, the profile function rapidly approaches the unitarity limit 
as x decreases when d n is large. This can be seen in Fig. EH where we have repeated the 
calculation of Fig. ^ this time using d„ = .8 fm. Note the large values of Th{b) at small b 
for d ~ .5 fm. Therefore, by choosing a smaller value for d n we are making a conservative 
estimate of the approach to the BBL. 



IV. ESTIMATING THE PROXIMITY OF THE TOTAL 7 *iV CROSS SECTION 
TO THE BBL 

To properly study the proximity of DIS to the unitarity limit, we must evaluate the degree 
to which the different hadronic components contribute to the 7*iV cross section, ctJ^t- T 
and L refer, respectively, to the transverse and longitudinal cross sections. In the color 
dipole formalism, the longitudinal and transverse cross sections can be factorized into the 
convolution of the perturbative light-cone wavefunction with a universal color dipole cross 
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FIG. 10: The hadronic configuration-nucleon profile function for different x values. Here we have 
used d n = .8 fm. Compare with Fig. ^ 



section, 



dz / d 2 d\i) L , T (z,d)\ 2 a tot (d,x) . 



(26) 



In this paper, Eq. (|26|) applies also to cases where at t{d, x) is the cross section for interactions 
of large size hadronic configurations with the nucleon. z is the quark momentum fraction, 
and ipL t T{ z ,d) is the longitudinal/transverse photon wave function calculated in QED. We 
have calculated a^ L ^ using a(d, x) with t-dependence determined in Sect. El Plots of F 2 
and Fl are shown in Fig. El and Fig. El The structure functions F 2 and Fl are defined as, 

Q 2 



F L (x,Q 2 
F 2 (x,Q 



4vr 2 



2 



-(<7 L + cr T ) 



(27) 



Note that the t-dependence is not needed for calculations of F 2 and F L . For more plots 
of the total cross section calculated within the QCD improved dipole model see Ref. 0. 
Figure 0] shows the distribution of the integrand in Eq. ()26|) over total hadronic sizes and 
demonstrates the suppression of large size hadronic configurations. 
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FIG. 11: The fraction of a^ L ^ with contributions coming from values of greater than the 
corresponding values listed on the x-axis. When Q 2 = 2 GeV 2 and x = 10 -4 , about 1/5 of <x 



7* AT 
t ot 



comes from hadronic components scattering with Th > 1/2. 



We would like to study the contribution of intervals of T(b) to the 7*iV cross section. 
However, the profile function for 7*iV scattering by itself is not useful because the photon 
wave function is not normalizable and because it depends on a e . m .. Thus, in order to look 
for the proximity to the BBL, we have plotted the fraction of the total 7*iV cross section 
due to different regions of the hadronic profile function. Plots with different values of x 
are shown in Fig. ^2 The y-axis denotes the fraction of the longitudinal (transverse) cross 
section with contributions from Th{b) greater than the corresponding value on the x-axis. 
Note that when Q 2 = 2 GeV 2 and x ~ 10 4 , about 1/5 of the longitudinal cross sections are 
due to dipole configurations corresponding to Th(b) > 1/2. 

At Q 2 = 20 GeV 2 , Fig. HU demonstrates the recovery of the leading twist behavior, 
especially for ai at x = .01, where less than one-tenth of o~l is due to hadronic configurations 
with Th(b) > 1/2. It is clear from Fig.^^that, for central impact parameters and low enough 
Q 2 , a significant portion of the total cross section is due to hadronic configuration-nucleon 
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interactions that are close to the unitarity limit. In high energy 7* A scattering, where 
the effects of the BBL are enhanced, we may be able to use DIS to probe the BBL. This 
possibility is discussed within the context of the QCD improved dipole model in Sect. IV II 



V. COMPARISON WITH RESULTS OF OTHER STUDIES 

The reasonableness of our model is demonstrated in Fig.lHlwhere the dipole model is seen 
to be consistent with recent HERA data for F2 at low values of Q 2 and x. Furthermore, in 
Fig. our model is seen to be consistent with preliminary results from HERA for Fi [26] . 
Other studies of the impact parameter picture of hadronic interactions with nucleons were 
done in where p production data was used to extract the t-dependence. The analysis 
in jsj used the S-matrix convention, S(b) = 1 — T(b), in place of the profile function. In 
Fig. El we have plotted our prediction of the S-matrix profile for central impact parameters 
along with earlier result from Munier et al. |2| . In their analysis, the authors were restricted 
to using p production data to model the t-dependence. Data for p production is limited 
to kinematics where —t < 0.6 GeV 2 (b > .3 fm) so the accuracy of their results is limited 
to moderate impact parameters. Further complication is introduced into their analysis by 
the need to model the p- meson wavefunction. Our model uses J /if) production data and is 
therefore valid at small values of b. As mentioned in Sect. HH production of J/tj) depends 
only on the two-gluon form factor, and can therefore be extended down to very small impact 
parameters. It is natural to compare our results for S(b) with the median value of the dipole 



size as it is was used in the calculation of the amplitude estimated in Ref. [27[. In our case, 
the value of d that corresponds to S(b) evaluated at Q 2 = 7 GeV 2 is about d ~ 0.32 fm 
With this assumption, Figure shows that our model has very good agreement with the 
results of Ref. |2j at moderate values of b while there is an expected deviation between the 

n 

two models for low values of b (see footnote |2j|). Since our model will have ~ t 4 behavior 

at small values of b, then even in the small b region, our model deviates from the results of 

Ref. I2J by no more than about twenty-five percent. Reference I2J used a simple exponential 

or power ansatz to interpolate to larger values of t as indicated in the figure. In Fig. it is 

n 

seen that the model used in |2j has a high degree of uncertainty at small values of b because 
of the necessity to guess the form of the function that interpolates to large t. In contrast, 
our model uses information about J / ip production and DIS to model the small b behavior. 
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FIG. 12: Comparison of the S-matrix calculated using Eq. IjlOj) with results obtained in 
bold line shows our result, while the dashed and dotted lines show results taken from Ref. 
p production with three different interpolations for the i-dependence. 
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using 



After preliminary results of our study were presented, there appeared an experimental 
analysis with improved data on inclusive cross sections and vector meson production at 
HERA by Henri Kowalski and Derek Teaney jy| . They carry out an analysis similar to that 
used in |2j. Therefore, it differs from our analysis in that it does not include information 
about large t behavior of the two-gluon form factor. Future improvements on the dipole 
picture should make comparisons with this data. 



VI. SCATTERING OFF A HEAVY NUCLEAR TARGET 

It is interesting to examine how the profile function T^s, b, d) changes when the free 
proton target is substituted by a heavy nuclear target such as the nucleus of 208 Pb. In the 
heavy nucleus case, the procedure for obtaining I\(s, b, d) differs from the one in the nucleon 
case and is outlined below. First, for dipoles of small transverse sizes, d < 0.2 fm, the inelastic 
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scattering cross section at a given impact parameter b is given by the perturbative QCD 
expression involving the impact parameter dependent nuclear gluon distribution, gA^x, Q 2 , b) 
(compare to Eq. (JTJ)) 



a^ CD (d, x, b) = ^-d 2 a s (Q 2 )xgA(x', Q 2 , b) , (28) 

where x' is given by Eq. (J23). The gluon distribution gA{x,Q 2 ,b), normalized such that 
J d 2 bgA(x, Q 2 ,b) = gA(x, Q 2 ), was evaluated in [gj using the theory of leading twist nuclear 
shadowing. The profile function I\(s,5, d) can be found from Eq. (|2*Tf (see see also 

2ReT h ( 8 , b, d) - \T h (s, b,d)\ 2 = a^ D (d, x, b) , (29) 

Ignoring the small imaginary part of r\(s, b, d), which is even smaller in the heavy nucleus 
case than in the free proton case because of the effect of nuclear shadowing, Eq. (|2H|) gives 



T h {s, b,d) = l- J l - at n Q ^ D (d, x, b) , (30) 



which is valid for d < 0.2 fm. Second, for dipoles of a larger size, do = 0.2 < d < d n = 0.65 
fm, the profile function is found by interpolating between the pQCD expression of Eq. (J30j) 
and the profile function calculated at d = d v j^J 

d 2 - d 2 

F h (s, b, d) = (T h (s, b, d n ) - T h (s, b, d )) + T h (s, b, d ) . (31) 

a n a 

The profile function Th(s,b,d n ) is calculated using the Glauber multiple scattering formal- 
ism [33] 

r h (s,b,d„) = 1 - e ^iv( S )T(&)/2 ^ (32) 

where a nN (s) is the energy-dependent pion-nucleon total scattering cross section, <t 7I -^y(s) — 
23.78(s/so) ' 08 mb, so = 200 GeV; T(b) is the nuclear optical density normalized such that 
J d 2 bT(b) = 1. A is the number of nucleons in the target. 

Third, for the dipoles with the size d > d n , the profile function is given by Eq. (|32j). where 
the pion-nucleon cross section is allowed to slowly grow as 

1 5 d 2 

a^ N (s, d) = a % N (s) d 2 + d 2, 2 - ( 33 ) 

The results for the profile function r^(s, b, d) for the nucleus of 208 Pb are presented in 
Fig. IIHlby the two solid curves (x = 10~ 4 and x = 10~ 5 ), dashed (x = 10~ 3 ) and dot-dashed 
(x = 0.01) curves. 

25 




-\ d: 


0.2 Cm 















5 10 
d=0.6 fm 




b, fm 



b, fm 




b, fm 



FIG. 13: The hadronic configuration-nucleus ( 208 Pb) profile function. The upper solid curves 
correspond to x = 1CT 5 , and the lower solid curves corresponds to x = 10~ 4 ; the dashed curves 
correspond to x = 10~ 3 ; the dot-dashed curves correspond to x = 10~ 2 . 

The profile function for the nuclear target shows some similarity with the profile function 
for the proton target. The main differences are that the BBL is approached over a larger 
range of impact parameters than in the case of a proton target. This is not surprising 
because of the larger thickness of the nuclear target. The plots in Fig. ITU show the fraction 
of the hadronic cross section due to large values of IV 

Large leading twist gluon shadowing tames the growth of the interaction of hadronic 
components of the photon with the nucleus so that the unitarity constraint is satisfied for 
x > 1CT 4 while the BBL may be reached for a large range of impact parameters. For smaller 
x, unitarity starts to break down at central impact parameters. For large d, unitarity is 
automatically satisfied since the Glauber model for large total cross sections leads to a 
that approaches unity. 
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FIG. 14: These plots are the analogue of those appearing in Fig. El They correspond to the profile 
functions for the nuclear target in Fig. ^3 In these plots, x corresponds to Q 2 = 2 GeV 2 . 

Finally, we have included plots in Fig. for the nuclear target showing the fraction of 
<t 7 * a due to large values of Th analogous to those in Fig. [TT] In Fig. HH1 we see that the 
BBL is approached for nearly all values of x at Q 2 = 2 GeV 2 and Q 2 = 20 GeV 2 . (Note the 
recovery of leading twist behavior for the longitudinal cross section at large x and Q 2 = 20 
GeV 2 .) Notice also that the fraction of a 7 *" 4 due to large values of Th for x = .01 is actually 
larger in some cases than for the case, x = .001. This effect can be explained qualitatively 
by inspection of Fig. ED For the case of the nuclear target, the main contributions to <j" 1 * A 
come from smaller values of d (d ~ .2 fm). The growth of the profile function with decreasing 
x at small d is slower for smaller values of x (x ~ .01) than for larger values. Thus, the tail 
of profile function at large impact parameter may become significant in these regions. 



VII. CONCLUSION 



A general, well-known feature of 7*iV and scattering is that the fraction of the 
interactions of hadronic components in the virtual photon wavefunction with the proton 
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FIG. 15: These plots are the analogue of those appearing in Fig 1111 The fraction of a^* A due to 
values of the hadronic profile function larger than Th is plotted versus T^. 

which take place at or near the BBL increases as x and Q 2 decrease as is exhibited explicitly 
in this paper in Figs. ^ El and ^2 We hope that one day we may be able to exploit the 
novel properties of interactions in the BBL to study a new phase of pQCD. In Ref. [34| the 
signatures of the BBL for DIS were discussed with the hopes that they would be seen in 
future experiments. It remained to be determined, however, in which kinematical regions one 
can expect to see black body behavior. Having now constructed a model of the amplitude 
for the interaction of the hadronic components of the virtual photon wavefunction, we are in 
a position to make rough estimates of the fraction of the hadronic interactions that exhibit 
the characteristic behavior of black body interactions. More precisely, since we know that 
the effects that we have ignored so far - inelastic diffraction and the real part of A^n-, - will 
tend to increase the proximity of the interactions to the unitarity limit, then we can place 
lower limits on the values of x and b where a significant fraction of the events will occur at 
or near the BBL. Our results show that, within available HERA kinematics, a significant 
fraction of the total DIS cross section is due to interactions of the hadronic components with 
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the proton that occur near the BBL. In particular, Fig. II II shows that at Q 2 ~ 2.0 GeV 2 and 
x < 10~ 4 , about 1/5 of the longitudinal cross section is due to values of Th(s,b,d) > 1/2. 
The agreement of our model with preliminary HERA data and with previous models helps to 
strengthen this conclusion. An improved model, with corrections for inelastic diffraction, will 
likely predict a more rapid approach to the BBL at small x and central impact paramters. 
The approach to the BBL as x and b decrease occurs much more rapidly for the case of a 
heavy nuclear target than for the case of a proton target. This can be seen by comparing 
Figs. 01 and ITU For example, at d — 0.4 fm and x ~ 10~ 4 , Fig. El shows that, for the proton, 
less than 1/2 of the total cross section is due to contributions from I\(s, b, d) > 0.5 whereas 
with a 208 Pb target, Fig. El shows that over seventy percent of the total cross section is 
due to contributions from T h (s, b, d) > 0.5. This suggests that nuclear targets are ideal for 
studies of the BBL as a phase of QCD as has been discussed before in Ref. js^J. Future 
work on this subject should incorporate inelastic effects. Also, a greater understanding of 
the large t behavior would lead to greater accuracy in the model. 
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